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movable, physically meaningful phases in quark mass matrices and (ii) which elements of 
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fermion generations; these provide further insight into the three-generation case of physical 
interest . 
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1 Introduction 

Understanding fermion masses and quark mixing remains one of the most important out- 
standing problems in particle physics. In an effort to gain insight into this problem, many 
studies of simple models of quark mass matrices have been carried out over the years. The 
phases in these mass matrices play an essential role in the Kobayashi-Maskawa (KM) mech- 
anism IQ for CP violation^. A given model is characterized by the number of parameters 
(amplitudes and phases) which specify the quark mass matrices. Thus, a very important 
problem is to determine, for any model, how many unremovable, and hence physically mean- 
ingful, phases occur in the quark mass matrices and which elements of these matrices can 
be made real by rephasings of quark fields. We recently presented a general solution to this 
problem ||^ . Here we give a detailed discussion of our results and proofs, and applications 
to currently viable models. To make the paper self-contained, we will review the results in 
0. In passing, it should be mentioned that we have also given a general solution to the 
analogous problem for lepton mass matrices |^0. The organization of this paper is as follows. 
In section 2 we discuss our theorem on the number of unremovable phases. In section 3 
we give the details of our complex rephasing invariants and results on which elements of 
the mass matrices can be rendered real. This section contains the proofs of certain theo- 
rems presented in Ref. 0. Section 4 contains some new results on complex invariants and 
unremovable phases for arbitrary Nq- In section 5 we apply our general results to models. 



^As will be discussed further below, there are, in general, more unremovable phases in the quark mass 
matrices than in the quark mixing matrix, so that some of the former phases contribute to CP-conserving 
quantities. We also note that our analysis does not assume that the KM mechanism is the only source of 
CP violation. 

^The situation in the leptonic sector is qualitatively more complicated than that in the quark sector 
because of the general presence of three types of fermion bilinears, Dirac, left-handed Majorana, and right- 
handed Majorana, each with its own gauge and rephasing properties. It is also different because one does 
not know the full set of fields which might be present: there might or might not exist electroweak singlet 
neutrinos. The most general case was treated in Ref. [0. The uncertainty in the field content of the neutrino 
sector produces the same uncertainty in a general analysis of phases for both the quark and lepton sectors 
considered together. An analysis of invariants and their constraints in a particular unified model was also 
given in Ref. j|]. 



2 Theorem on the Number of Physical Phases 

The quark mass terms are taken to arise from interactions which are invariant under the 
standard model gauge group Gsm = SU{3) x SU{2) x f/(l), via the spontaneous symmetry 
breaking of Gsm- In the standard model and its supersymmetric extensions, the resul- 
tant mass terms appear at the electroweak level via (renormalizable, dimension-4) Yukawa 
couplings^'''. It follows that these mass terms can be written in terms of the Gsm quark 
fields as 

-C^=Y, [{Q,L)lM%KkR + (Q,L)2M,?4ij] + h.c. (2.1) 

where j and k are generation labels; Nq = 3 is the number of generations of standard-model 
fermions; Qjl is an SU(2) doublet, with Qil = {d)L, Q2L = {^s)l, Qzl = (Dl] the subscript 
a on {QjLla is the SU(2) index; and the SU(2)-singlet right-handed quark fields are denoted 
as UkR with uiR = ur, U2r = cr, u^r = Ir, diR = (Ir, ^r = sr, d^R = Br. As indicated, 
we will concentrate on the physical case of A*",:; = 3 generations of standard-model fermions 
(with associated light neutrinos). However, for the sake of generality, and because it provides 
further insight into the physical Nq = 3 case, we will give a number of results for arbitrary 
Ng- M*^"-* and M^^^ are the mass matrices in the up and down sectors, whose diagonalization 
yields the mass eigenstates Ujm and djm- Because of the assumed origin of the quark mass 

■^We do not consider models in which fermion masses arise via multifermion operators, which are not 
perturbatively renormahzable. Note that at a nonperturbative level, lattice studies H show that a lattice 
theory with a specific multifermion action and no scalar fields may yield the same continuum limit as a 
theory with a Yukawa interaction. 

■^In a number of interesting models, these Yukawa couplings are viewed as originating at a higher mass 
scale, such as that of a hypothetical supersymmetric grand unified theory (SUSY GUT) or some other 
(supersymmetric) theory resulting from the E << (a')^^" limit of a string theory. It is possible that 
interactions which appear to be Yukawa couplings at a given mass scale, could be effective, in the sense that 
some of the elements of the associated Yukawa matrices could actually arise from higher-dimension operators 
at a yet higher mass scale. An example of a higher-dimension operator which might plausibly occur near 
the scale of quantum gravity is a dimension-5 operator of the generic form (c/Mp)0i02Qi9j,, where (f)i and 
(1)2 are appropriate Higgs fields, c denotes a dimensionless coefficient, and Mp is the (reduced) Planck mass, 
Mp = ■\/?ic/(87rGjv). The plausibility of such higher-dimension operators can be inferred either from the 
nonrenormalizability of supergravity or as a consequence of the E << {a')~^'^ limit of a string theory. At 
lower mass scales this operator could produce effective Yukawa terms contributing to (2.1), where, say, (f>i 



contains Hi, i = 1 or 2 as defined in (2.3), (2.4), and the factor c < 02 > /Mp enters into the effective 
Yukawa or mass matrix (where < (f)2 > might have a value somewhat smaller than Alp but much larger than 
the electroweak scale). 



matrices from Yukawa interactions (i.e. interactions which appear as dimension-4 at a given 
mass scale, but which could include the contributions of higher-dimension operators at higher 
mass scales) , it is convenient to write these mass matrices in terms of dimensionless Yukawa 
matrices according to 

F(/) = 2'^%/M^f^ , f = u,d (2.2) 

where Vu and Vd are real quantities with the dimensions of mass, and represent Higgs vacuum 
expectation values in the standard model and its supersymmetic extensions. For example, 
in the minimal supersymmetric standard model (MSSM), 

< ffi >= ( ,\, ] (2.3) 



and 



<H2>=( 2"'q'^" ) (2.4) 



where Hi and H2 are the Y = 1 and Y = —1 Higgs fields. With one definition for the angle 
P in the MSSM, v^ = f sin/?, Vd = f cos/?, where v = 2^^^^Gp = 246 GeV (many authors 
use a different convention according to which t>„ = t>cos/?, Vd = f sin/3). Because f„ and Vd 
are real, the phase properties of the mass matrices M^-^^ are the same as those of the Yukawa 
matrices Y^^^ and we shall deal with them interchangeably.^ 

To count the number of unremovable, and hence physically meaningful, phases in the 
quark mass matrices, we rephase the fermion fields in ( |2.1| ) so as to remove all possible 
phases in these matrices. 

A general theorem on this counting problem was given before 0; we review it here. 
As was evident in Ref. 0, this theorem applies for arbitrary A^^; although the physically 
interesting case is Nq = 3, we indicate the general A^c-dependence here. One can perform 
the rephasings defined by 

Q^L = e-^'^^Q'^L (2.5) 

^In specific models in which the Yukawa couphngs arise at a high mass scale, e.g. in a hypothetical grand 
unified theory or an effective pointlike field theory arising from the E << (a')~^" limit of a string theory, 
and in which below this scale the gauge group reduces to Gsm without further intermediate structure, it 
is the y'-'i^ which one evolves down to the electroweak scale using the appropriate renormalization group 
equations. Note that M^^^ and Y^^^ as defined above may subsume the effects of several different types of 
Yukawa couplings of fermion and Higgs representations in the theory at the high mass scale, e.g. a GUT. 



u,R = e'^'^'v!^^ (2.6) 

d,K = ^^^d:^K (2.7) 

for 1 < j < A^G- (The minus sign in ( ^.5[ ) is included for technical convenience to avoid 
minus signs in later formulas.) In passing, we note that since the mass terms in eq. (|2.1| ) are 
assumed to have arisen from G^jif-invariant Yukawa interactions, the Yukawa matrices, as 
coefficients of GsAf-invariant operators, are obviously left invariant by Gsm gauge transfor- 
mations. Also, note that since at the level of ( |2.1| ) which defines the quark mass matrices, 
no explicit Higgs fields appear, we do not apply the rephasings of these Higgs fields. In any 
case, in particular models, such as the MSSM, the original Higgs doublets ii\ and H^ have al- 
ready been rephased so as to yield real vacuum expectation values v^ and f„. In terms of the 
primed (rephased) fermion fields, the mass matrices, or equivalently the Yukawa matrices, 
have elements 

yU)' ^ S^M'^K) (2.8) 

for f = u,d. Thus, if Y^^'^ has Nf nonzero, and, in general, complex elements, then the Nf 
equations for making these elements real are 

a,+/3i^^ = -argiYJP)+v%^7r (2.9) 

for f = u,d, where the set {jk} runs over each of these nonzero elements, and r]A = or 
1. The rijk term allows for the possibility of making the rephased element real and negative 
rather than positive; this will not affect the counting of unremovable phases. Let us define 
the vector of fermion field phases 

^ = ({«.}, {/??^},{/?f^}r (2.10) 

of dimension 3Ng, where {aj} = {ai, ...,aNG}y {Pi } — {Pi \ ■■■yPNc} ^°^ f = '^yd, and 

w = ({- arg(y/:)) + r/g^vr}, {- arg(llt)) + r/^^) tt})^ (2.11) 

of dimension equal to the number of rephasing equations Nf.q = N^ + N^, where the in- 



dices jk and mn take appropriate values]^ For explicit calculations, the ordering of the 
elements of the vector w is given by the sequence of nonzero elements Y J^ for jk = 
{11,12,13,21,22,23,31,32,33} followed by the nonzero elements YJ;^^ with mn = ranging 
over the same set, {11, 12, 13, 21, 22, 23, 31, 32, 33}. We can then write (|2.9| ) for f = u,d as 

Tv = w (2.13) 

which defines the N(,g-Tow by 3A'"G-column matrix T. Since a nonzero element of Y^-l'\ 
f = u,d is, in general complex (because no symmetry requires it to be real), there is an 
obvious general relation between the number of rephasing equations N^g, its maximal value, 
{Neq)max = '^^Gi ^"^^ ^^^ numbcr of zero elements, iV^, in the Yukawa matrices: 

iVe, = 2iV2 - iV, (2.14) 

In studying specific models, it is also helpful to introduce the notation N^j for the number of 
zeroes in each of the Yukawa matrices Y'^^\ for f = u and f = d. Obviously, A^^ = N^^u+N^^d- 
If an element of Y^-f\ f = u, d is zero at some energy scale E, this will not necessarily be 
true at a different energy scale E', so that in specifying the numbers N^^u and N^^d, as well as 
in specifying the forms of the Yukawa matrices F*^"^ and Y^'^^ themselves, one must indicate 
the energy scale involved. If at this energy scale, the gauge group and the structure of 

the Yukawa couplings guarantee that the Y^-f\ f = u,d, are (complex) symmetric, then, of 

( f) (f) 

course, Y^^ = is equivalent to Y^j = 0; however, it will be convenient to count both in 

our definition of N^j and N^^. 

Since T is an Neq-row by 3A^G-column matrix, clearly 



ran 



k{T) < min{iV,„ 3Ng} (2.15) 



6t 



In eq. ( 2.11 ) we chan ge notation from Ref. y], where an overall minus sign was absorbed in T, so that 



the analogue of eq. ( 2.11 ) was 

w = ({arg(F(")) + ,yj:V}, {arg(r('^)) + ^f^n}f (2.12) 

(recall that since the arguments are defined mod 27r and 77^-^ = or 1, it follows that rjji. tt — —rj/k t^ (mod 
27r).) The present definition avoids minus signs in our explicit listings of T matrices below (e.g. eq. ( ^.1.3 ) 
etc.). 



For the case (realized in all experimentally viable models which we have studied) where 
Neq > 3Ng, one has a more restrictive upper bound on rank{T): 

rank{T) < 3Ng - 1, {for N^g > 3Ng) (2.16) 

This is proved by ruling out the only other possibility, i.e. rank{T) = 3Ng- The reason that 
rank{T) cannot have its apparently maximal value is that one overall rephasing has no effect 
on the Yukawa interaction, namely the U(l) generated by (|2.5|) -(^l7|) with — «« = l3j = Pf, 
for all i,j,k. For the relevant case Nq = 3, for all of the (realistic) models that we have 
studied, the upper bound ( |2.16| ) is saturated, i.e. rank{T) = 8. This is not a general result; 
we will exhibit a toy model (model 7 in section 5) where rank{T) = 7, although that model 
is a degenerate case. 

The main theorem P| is as follows: 
Theorem: 
The number of unremovable phases Np in Y^^^ and Y^'^^ is 

Np = Neq-rank{T) (2.17) 

(which holds for arbitrary No)- 

Proof: 

First, denote rank{T) = r^. Then one can select and delete Nf,g — r^ rows from the matrix 

T without reducing the rank of the resultant matrix. This deletion means that we do not 

attempt to remove the phases from the corresponding elements of the Y^-f\ f = u,d. We 

perform this reduction. For the remaining r^ equations, we move a subset of SNc — r^ phases 

in V to the right-hand side of ( |2.9| ), thus including them in a redefined w. This yields a set 

of Tj, linear equations in r^ unknown phases, denoted v. We write this as Tv = w. Since by 

construction rank{T) = r^, it follows that T is invertible, so that one can now solve for the 

r^ fermion rephasings in v which render r^ of the N^g complex elements real: v = T~^w. 

Hence there are Neg — r^ remaining phases in the Y^-l'\ as claimed. □. 

The operations involved in this proof and the construction of T will be explicitly illustrated 

for model 1 in section 5. 
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Some comments are in order. First, we note that in order for these phases to be physically 
meaningful, it is, of course, necessary that they be unchanged under the full set of rephasings 
of fermion fields. In fact, as we will show, there is a one-to-one correspondence between each 
such unremovable phase and an independent phase of a certain product of elements of mass 
matrices which is invariant under fermion field rephasings. (In special cases a model may 
have an unremovable, invariant phase which, because of a particular symmetry, is zero or tt.) 
Second, as is clear from our proof, in general, the result ( p.l7| ) does not depend on whether 
or not Y I = YH initially. Hence making Y^^^ (complex) symmetric for either f = u or 
f = d does not, in general, result in any reduction in Np, since in terms of the Gsm theory, 

nothing guarantees this symmetry, and the rephasing is carried out in terms of Gsm fields. 

( f) 
Third, if one of the unremovable phases is put in a given off-diagonal Y ): , one may wish to 

modify the kj'th equation to read 

Ok + /?f = - arg(Fif ) - arg(r//)) (2.18) 

For example, in a model where \Y-l \ = \YM'\, this would yield Y I = Y"^ for this pair jk. 



The modification in ( 2.18| ) has no effect on the counting of phases. 



Using our theorem, we may calculate the maximal number of unremovable phases Np 
in the Yukawa matrices as a function of the number of fermion generations, Nq. For this 
purpose, we consider the case in which these matrices have all nonzero arbitrary (and, hence 
in general, complex) elements. Then N^^g takes on its maximal value, 

(iVeJmax = 2N^ (2.19) 

Since the elements of the Y^^^ are arbitrary, it follows that they satisfy no special relations, 
so that the T matrix has maximal rank. By our result (|2.16|) , this is rank{T)fnax = 3A^g ~ 1- 
Using our theorem ( p^.lT] ), we thus find that the maximal number of unremovable phases is 

{NpUa. = {Ng - 1){2Ng - 1) (2.20) 

That this case yields the maximal number of unremovable phases is clear, since we can 
obtain any other form by making some entries of the Yukawa matrices zero, and this can 



only decrease the number of unremovable phases. For Nq = 1,2 and 3, {Np)max takes the 
values 0, 3, and 10, respectively. In current phenomenological models of quark mass matrices, 
one tries to make as many entries in the Yukawa matrices as possible zero to minimize the 
number of parameters and increase predictiveness, so that in the physical Nq = 3 case, we 
will find that Np is typically just 2 or 3. 

3 Rephasing Invariants and Theorems on Locations of 
Phases 

A fundamental question concerns which elements of y-"-* and F*^'^-' can be made real by 
fermion rephasings. This is connected with the issue of which rows are to be removed from 
T to obtain T, i.e. which nonzero elements of the F*^"^ and y'-'^-' are left complex. In Ref. 
several theorems were given which provide a general answer to these questions. We 
review these theorems here and give the details of the proofs. As in the previous section, 
for generality, we shall indicate which results apply for an arbitrary number of fermion 
generations, Nq- In the next section we shall also give a number of new results for arbitrary 
Ng- The general method is to construct all independent complex products of elements 
of the Y^f\ f = u,d, having the property that these products are invariant under the 
rephasings ( |2.5| )-( P?^ ). Here we use the term "invariant products" to mean invariant under 
the rephasings ( |2.5| )-( P77| ). These products must involve an even number of such elements, 
since for each index on a Y, there must be a corresponding index on an Y* in order to form a 
rephasing-invariant quantity. Since in general, by construction, these invariant products are 
complex, i.e. have arguments 7^ 0, vr, each one implies a constraint which is that the set of 2n 
elements which comprise it cannot be made simultaneously real by the rephasings ( |2.5| )-( pr7D 
of the quark fields. Here, by "complex invariant products" , we mean products which are, in 
general, complex for arbitrary Y^'l'^; in special cases, a symmetry of a given model may render 
some of these complex invariants real. It is easily seen that the only rephasing-invariant 
products of two elements of the y*--^-* must be of the form 1^/ P, i.e., there are no complex 
invariant products of order 2. We define an irreducible complex invariant to be one which 



cannot be factorized purely into products of lower-order complex invariants. Given that there 
are no quadratic complex invariants, it follows that all complex invariants of order 4 and 6 are 
irreducible. As we will prove below, for the physical case of Nq = 3 generations of fermions, 
the phase constraints are totally determined by the complex invariants of order 4 and 6, so 
in the analysis of phenomenological models, one automatically deals only with irreducible 
complex invariants.]] We define a set of independent (irreducible) complex invariants to be a 
set of complex (irreducible) invariants with the property that no invariant in the set is equal 
to (i) the complex conjugate of another invariant in the set or (ii) another element in the set 
with its indices permuted. This does not imply that the arguments of a set of independent 
(irreducible) complex invariants are linearly independent. Indeed, we will show examples 
below of sets of independent complex invariants of the same order whose arguments are 
linearly dependent, and also examples of higher-order complex invariants whose arguments 
can be expressed as linear combinations of arguments of lower-order complex invariants. 
We define Nia to denote the number of linearly independent arguments {—>■ ia) among the 
Ninv independent complex invariants. Further, we define Ninv to be the total number of 
independent (irreducible) complex invariants in a given model. It is useful to define Ninv,2n 
to denote the number of independent complex invariants of order 2n. 

We thus construct a set of rephasing-invariant products depending on the up and down 
quark sectors individually: 

n 
pif) ^ TT y (/) yif)* (o i N 

a=l 

where f = u,d, and ai is an element of the permutation group Sn- (This set of invariants is 
clearly the same as the sets defined by JlLi ^/il^'/i^Cfc.) and Ua^i YlSjiuhanika)^ ^^^^re 
c"_R ^ Sn', with no loss of generality, we write it in the form of eq. ( p.l|) .) The first theorem 
is: The products in { \i.l\ j are invariant under the rephasings of fermion fields ( ^.5D -( PT7D . To 

^Reducible complex invariants only occur in the hypothetical case Nq > 4 and then first occur at order 
8, as products of different complex quartic invariants. Since reducible invariants do not yield any new phase 
constraints, it suffices to consider only irreducible complex invariants. Hence in all of our discussion below, 
in particular, the section on arbitrary Nq, we shall take "complex invariants" to mean "irreducible complex 
invariants" . 



prove this, we observe that from eq. (|2.8|) , it follows that under these rephasings, 

a=l a=l 

n 

V / 11 Jaka-^ CTLija)ka 

a=l 
n 

= n ^jaka^aL{ja)ka (3-2) 

a=l 

where Zlb=i(«ii, — <^o-l(j6)) — because the permutations of Sn are an automorphism of the 
set (1,2, ...,n). Hence, the products -P2n;jiA:i,...jnfen;o-L ^^^ invariant under rephasings of the 
quark fields, as claimed. □. 

Secondly, we construct a set of invariants connecting the up and down quark sectors: 



s t s t 

^2n;{j},{k},{m};aL,au;cTa Vll jakaf\l.L ^ js+b "itAli ^ aL{jc)cTu{kc) J \1.L o-i(js+e)o-d{me)^ 

a=l 6=1 c=l e=l 



where s, t > 1, s + 1 = n, ctl G S'„, cXu G S's, and 0"t G St- The proof that these are invariant 
proceeds in a manner similar to that just given above; hence we do not give it explicitly. 
At quartic order, 2n = 4, cr^ G 5*2 in eq. (|3.1|) . We use the standard notation (Jj) = 1 and 



(2^) = (12) = T (transposition) for the elements of 5*2. For a^ = t, we obtain the complex 
invariants 

p{f) ^ pU) 

jiki,J2k2 ~ 4,;jiki,J2k2;<rL=T 



_ vif) vU) \rif)*T^{.f)* /q A\ 

~ ^jlki^J2k2^J2ki^jik2 y^-^) 

for f = u,d. (For the other choice, a, = 1, the resultant invariant product, Pi£,,,,,,.,=i = 
l^/iiiPl^jlfel^ is real and hence yields no constraint on phases.) At this quartic order, there 
is only one Q-type complex invariant; this has s = t = 1, o"l = r, so that in the general 
notation of (PT^), it is Q^^.'^-^j^y^^^^y^^^^y.^^^^^^^^^^^^^-^. For brevity of notation, we denote it 
simply as 



^jiki,J2mi — V4;{jij2},{A:i},{mi};o-L=-r,tTu=l 



y{u) y{d) yH*y(d)* /q CT \ 

jlki J2mi-'- J2ki ^ jimi K'^-^J 
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(As with the quartic P invariant, the other choice for the permutation, 0"^ = 1, yields the 
product I^X^lmJ^' which is real.) Note the index sjTiimetry 

pif) ^ pif) (o c) 

jiki,J2k2 J2k2,jiki \'~'-^J 

and the complex conjugation relations 

pif) — p(f)* I'o -tX 

jikiJ2k2 jik2,J2ki \^-') 

and 

At order 2?2 = 6, we find one independent P-type complex invariant for each quark sector 
f = u,d, and two independent Q-type complex invariants; these were given in Ref. 0. Here 
we present the details of the analysis. For a given /, we construct the invariants correspond- 
ing to each of the six permutations ai G S3. In standard notation, these permutations are 
listed as 

S3 = {(III) = 1, (III) = (12), (III) = (23), (If,) ^ (13), (III) ^ (231), (III) ^ (312)} (3.9) 

Clearly for ai = 1, ^6y!fcij2fc2,i3fc3;^L=i = I ^/i fc 1 ^iS ^i£ I ^ ^^^^^ ^^ ^^^^ ^^^ ^^^^^ implies no 
constraint on phases. For the three transpositions ai = (12), (23), (13), the resultant product 

( f) 
reduces to a real factor times a quartic P invariant; for example, -Pg /ifci 72^2 j3k3-aL=(i2) ~ 

l^hk^l'^ ■^i!jiki,j2k2^ ^^^ ^'^ forth for the others. The last two elements of 5*3 yield genuine 6'th 
order complex invariants: 

pif) ^ v'(^) v^f^ v^f^ Y^f>v^f>v'^f> rq 1 n^ 

-^6;iifcij2A;2j3fc3;o-L = (231) ^ jiki-^ J2k2-^ jaks-^ J2ki ^ j:jk2 ^ jik^ \o.lk}) 

and 

pif) ^ v'(^) v'^f') v^f^ Y^^f^v^f^v'^f^* I"? in 

^6;nki,J2k2,J3k3;aL = i312) ^ jiki^ J2k2^ jsks^ jiki ^ jik2 ^ J2ks V'^-^^) 

However, only one of these is independent, since they are related according to 



pif) _ pif) (n 1 2^ 

-^6;iifcij2fc2j3fc3;(TL=(231) -^6;i3fcijifc2,i2fc3;fTL = (312) W'-L^J 
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This explicit calculation proves our assertion above, that there is one (independent) 6'th 
order P-type complex invariant. We may take this to be the one in eq. (|3.11| ), and we 
denote it in a simple notation as 

pif) ^ pif) (o io\ 

^ jiki,J2k2,J3k'j, — ^ G-jikuJ2k2,J3k3;'TL={312) yo.io) 

To find the 6'th order complex Q-type invariants, we consider the case {s,t) = (2,1) 
in eq. (|3l3|); the other case (s,t) = (1,2) can then be obtained easily by the replacement 
y{w) ^^ y((i) [t^ f^]^Q resultant products. Since ai G S3, cr„ G 5*2, and ad = 1, it follows that 
there are 3! x 2! = 12 invariants of this type. Each of these can be labelled by cl ® (T„. 
Clearly the 10 1 and (12) (12) products are real. Further, we find that the products 
1 C?) (12), and (12) ® (12) give real factors times quartic P-type invariants, while (23) ® 1, 
(13) C?) 1, (231) (8) (12), and (312) ® (12) give real factors times quartic Q-type invariants. The 
remaining four choices for the permutations give genuine complex 6'th order invariants. To 
discuss these, we shall introduce a convenient notation, 

^e;jiki,J2k2,J3mi;aL,cTu = ^6;{ii,i2,J3},{fci,fc2},{mi};fTL,(Ti,{(7d=l) I'J-l^j 

The permutations (23) ® (12) and (312) ® 1 yield the same product, 

^(2,1) _ -i^H -iaH -iaW ■^(M)*y-(u)*-^(d)* /q 1 r\ 

^e;jiki,J2k2,J3mi;(7L={23),an = (l2) ^ jiki^ J2k2-^ J3mi^ jik2 ^ J3ki ^ J2nii yo.LO) 

while the choices (13) ® (12) and (231) ® 1 yield the same product, 

^(2,1) _ -i^(«) -i^C") -i^W v^{«)*-i^(«)*-iaW* ("?1fi'l 

^e;jiki,J2k2,J3mi;aL={13),au = il2) ^ jiki^ J2k2^ J3mi^ jgki ^ J2ki ^ jimi yo.lV) 



However, (|3.15|) and (|3.16| ) are not independent, since 



^6;iiA;ij2A:2j3mi;CTL=(13),(T„ = (12) Q6;J3ki,jik2,J2mi;aL = {23),au={12) (3-17) 

Hence, as stated above, there are only two independent complex 6'th order invariants, which 
may be taken to be (|3.15|) and its mirror pair under Y^"^^ <-> Y^'^\ For brevity of notation, 
we shall denote these simply as 

nifff) = yif) yif) yif) y(/)*y(/)*y(/')* /o 1 o\ 

^jiki,J2k2,J3mi jiki^J2k2 J3mi-^jik2 jski J2ini {O.LO) 
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where (fff) = (uud) or (ddu). Note the relations 

pif) ^ pif)* /o iq\ 

Q\jiki,J2k2,j'ik'i J3k2,J2ki,jiks \^u.j-^; 

and 

^6;jiki,J2k2,J3mi ^6;jik2,J3ki,J2mi yj.i^Kj) 

So far, all of the results in this section hold for arbitrary A^^. For the physical case Nq = 3, 
the order 4 and order 6 invariants given above form a complete set, since, as was noted in 
Ref. [^ (for Nq = 3) there are are no new constraints from any invariant of order > 8 (see 
below). We shall sometimes refer to P and Q invariants of order 4 and order 6 generically 
via the respective symbols P4, Q^, Pq, and Qq. 

We next review our theorems p[ and present the proofs. For a given model, construct 
the maximal set of Nia independent complex invariants of lowest order (s), whose arguments 
(phases) are linearly independent. Then we have the 
Theorem: 

(a) each of these invariants implies a constraint that the elements contained within it cannot, 
in general, all be made simultaneously real; (b) this constitutes the complete set of constraints 
on which elements of Y^'^^ and Y^'^^ can be made simultaneously real; and hence (c) Np = N^a- 
These results hold for arbitrary Nq- Further, (d) for the physical case Nq = 3, the complex 
invariants of order 4 and 6 form a complete set, i.e., the argument of any complex invariant 
of order > 8 can be expressed in terms of the arguments of complex invariants of orders 4 and 
6, and hence these higher order complex invariants do not yield any new phase constraints. 
Hence also, for the physical case Nq = 3, the total number of independent (irreducible) 
complex invariants is given by 

A^m. = Ar,„,,4 + Ar,„,,6, for No = 3 (3.21) 

(A discussion of complex invariants in the case of arbitrary A*",;; is given in section ^) 
Proof: The proof of (a) is essentially obvious: since by construction, the argument of each 
complex invariant is invariant under the rephasings of the fermion fields ( |2.5| )-(2r7) and since 



it is not, in general, equal to or vr, the corresponding invariant cannot in general be made 
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real. It follows, a fortiori, that it is not, in general, possible to make the Yukawa matrix 
elements which comprise this complex invariant product all simultaneously real. Given that 
this is a maximal set, the assertions in parts (b) and (c) of the theorem follow immediately. 

The proof of (d) is more involved. Let X2n be some non-zero invariant of type P or Q. 

( f) 
Further, let Y^j , f = u or d, he one of the elements in X2n- Then it is clear that in order for 

(f)* ( f)* 
X2„ to be an invariant, it must also contain the product of Y-)^ Y^^ for some fc, / = 1, 2, 3, 

where /' may be equal to, or different from, /. Hence we can write 

-^2n — l-lj llj i-ik ^ \6.ll) 

Now we consider the following two cases, (i) For some values of k and /, the Yukawa matrix 
element Yi): is non-zero. Then 

arg(X2„) = s.TgiYi^'>Y,P X^n) = 
^^siY^* (Yry!/^Y!i'>Y}/>) X') = arg{X, X,^^,) (3.23) 



where 



and 



Hence, 



X. = YrY^^Yii'>Y^> (3.24) 



X2n-2 — Yi^X (3.25) 



arg(X2„) = arg(X4) + arg(X2„_2) (3.26) 

Observe that eq. (|3.24|) subsumes both the case where X4 is of P type and the case where 



it is of Q type. Given the definition ( pJ.22D , it is easily seen that ( |3.25| ) is an invariant. 



Therefore, in this case the argument of the order 2n invariant X2„ can be expressed as a 
sum of the arguments of an order 4 invariant X4 and an order (2n — 2) invariant X2„_2, so 
that arg(X2„) does not yield a new constraint. By induction, every invariant of order greater 
than 4 which satisfies the conditions of (i) is reducible in this manner. 

It remains to consider the other case, that (ii) for all values of i, j, k and /, Y";^ =0. In 
this case, the invariant X^n cannot be reduced by means of the above procedure. Then it is 
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( f') 

also true that since Yi): vanishes, it cannot be one of the elements in the X2„, because the 

latter product is nonzero by assumption. Consequently, the invariant X2„, which contains the 

( f) ( f)* (f)* / (f) 

product Y^j Yij Y^)^ X , cannot contain the element 1^^ = at the same time. Because 

each invariant contains either no elements, or at least two elements, from each column of the 

Yukawa mass matrices (in the latter case, denote the two as Y^p and Y^p*), it follows that 

the X2n invariant in this case must be constructed out of blocks B^{^1 = (Y^pY^p*) such that 

the pairs of the first indices, (p, g), of the two such blocks do not coincide: 

X2n = n^i{ir> where {pi,qi) ^ {Pj,qj) if i^j (3.27) 

i 

Since p and q can take only 3 values each, there are (2) = 3 possibilities to assign these 
indices (where (^) = m\/{n\{m — n)\) denotes the binomial coefficient). This means that the 
order of an invariant which is not decomposable by these means cannot be greater than 3 
times the order of B, i.e., 6. We can exhibit examples (see section 5) for which this bound is 
saturated, i.e. the order of the lowest-order nonvanishing complex invariant is 6. This thus 
proves (d). □ 

For most models the lowest-order nonvanishing complex invariants are of order 4, but, as 
mentioned, we will in section 5 exhibit some exceptions for which the lowest-order nonvan- 
ishing complex invariants are of order 6. 

Corollary 1: 
If in a given model there are as at least as many independent quartic complex invariants 
with independent arguments as there are unremovable phases, then the arguments of not 
only order 8 and higher, but also of 6'th order complex invariants are expressible in terms 
of those of the quartic invariants and hence yield no new phase constraints. 
Proof: 

This is obvious, since the premise of the corollary means that since all of the Np = Nia 
unremovable phases have already been accounted for by quartic complex invariants, there 
cannot be any further phase constraints from any higher order complex invariants, and hence 
the arguments of all such higher order complex invariants must be expressible in terms of 
those of the Nia quartic complex invariants with independent arguments. □ 
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For most models, the lowest-order complex invariants are of order 4, but, as mentioned, we 
will exhibit two exceptions in section 5, where the only complex invariants are of order 6. 

Corollary 2: 
If in a given model all of the elements of the Y^^\ f = u,d are nonzero arbitrary (complex) 
quantities, then the arguments of the order 6, as well as the higher order, complex invariants 
are expressible in terms of those of the quartic invariants. 
Proof 

This is clear, since given any order 6 complex invariant, we can always apply the reduction 
method of eqs. (|3.23| )- (|3.26|) to express its argument in terms of that of the arguments of 



two quartic complex invariants. □ 

(Recall that in this case such a model would also have the maximal number of unremovable 

phases, (|]20D.) 

We have developed a graphical representation for the complex invariants which we have 
found useful. This is described next. The Yukawa matrix Y^^^ is a Nq x Nq matrix whose 
rows and columns are labelled respectively by the first and second indices, j and k on the 
elements K,/ . From the defining equation, ( p.4|) , we can represent the quartic invariant 



Pj!Lj2k2 = ^hkJjSjjS*Yj!k2 as follows. To the factor YJflyl^^* we associate an oriented 
line segment (which can be visualized as a vertical arrow t or | connecting the two indicated 
elements at the intersections of the rows ji and J2 with the column ki. Let the head of the 
arrow or equivalently the of the dipole be placed on the unconjugated element, and the tail 
of the arrow or equivalently the ® of the dipole be placed on the conjugated element. For the 
figures, it is actually convenient to use a notation involving vertical dipoles 1=1 and 1=1, 
where, as indicated, the of the dipole corresponds to the head of the arrow (unconjugated 

element), while the of the dipole corresponds to the tail of the arrow (conjugated element). 

(f) (f)* 
To the factor Yj2k^hk2 ^^ associate a similar oriented vertical line segment connecting the 

elements at the intersections of the J2 and ji rows with the column ^2. Note that the head of 

this arrow points in the opposite direction from the first arrow, or equivalently, the dipole is 

oppositely oriented. These two opposite arrows link the same rows; consequently, they form 

a rectangle. Hence, we may associate with each complex quartic P^-^-* invariant a (possibly 
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square) rectangle in the matrix Y^^\ As an example, our graphical representation of the 
quartic invariant P2233 is illustrated in Fig. 1(a). 

Similarly, we can represent the quartic invariant Qjiki,j2mi = ^/ifci V2mi ^2^1 iimi ^^ 
follows. Imagine the Y^"^' matrix to lie next to the Y^'^' matrix on a plane. To the factor 
^jifci jlfei *^ associate an oriented line segment linking the elements at the intersections of 
the ji and J2 rows with the kith column in the Y^^^ matrix. As before, the head of the 
arrow is assigned to the unconjugated element. To the other factor K,;,m^Y^^^„* we associate 
an arrow linking the elements at the intersections of the J2 and ji rows with the mi column 
in the F*^'^^ matrix. As with P, the two arrows making up the Q point in opposite directions. 
We may thus associate with each complex quartic Q invariant a rectangle linking the arrow 
in the Y^"^' matrix with the oppositely oriented arrow in the Y^'^' matrix. As examples, 
the quartic invariants (5i2,22 and Q23,32 are illustrated in Fig. l(b,c). The operation of 
complex conjugation of a given complex quartic P ot Q invariant reverses the directions of 
the arrows in each of the line segments. Since the complex conjugate invariants give the same 
phase constraints, the corresponding graphs with reversed arrows may be considered to be 
equivalent to the graphs already considered. Hence for purposes of enumerating invariants, 
one may suppress the degree of freedom associated with an overall reversal of all arrows 
(given that the pairs of columns forming each graph always have oppositely directed arrows), 
and thus deal with just the rectangles. This is symbolically indicated in Fig. 2 using the 
equivalent dipole notation. 

Higher order invariants may be represented in a similar manner. A 6'th order P complex 
invariant in the sector / (= -u or d) is represented as three vertical arrows in three columns of 
y(-^) with the property that the head of each is matched by the tail of some other arrow. A 
6'th order Q(""'^) complex invariant is represented by two oppositely directed vertical arrows 
in y*^") together with one vertical arrow in Y^'^\ again with the property that the head of 
each arrow is matched by the tail of some other arrow. Finally, a 6'th order Q('^'^") invariant 
is represented as for Q(""'^) but with the interchange of Y^^^ and Y^'^\ As examples, we show 
Q32S12 and Pii,23,32 in Figs. Figs. 3(a) and 3(b). 

Since the full set of Ninv independent complex invariants will have arguments which are 
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not, in general, linearly independent, it follows that 

Ninv > Nia (3.28) 

It may also happen that, e.g. for order 2n = 4, the number of independent quartic complex 
invariants, Ninv,4. is greater than N^a- For each complex invariant of a given order, X, 

arg(X)= Y: E4Sarg(F//)) (3.29) 

f=u,d j,k 

where the sum is over the Neq complex elements of Y^'^' and Y^'^'. Including all orders of 
invariants, these equations can be written as 

^ = Zw (3.30) 

where ^ is the A^j„t;- dimensional vector 

e=(arg(Xi),...,arg(X;v_))^ (3.31) 

Z is an A^j„^-row by A^gg-column matrix, and it is implicitly understood here that all of the 
r]-^ in eq. ( p.ll| ) for w are taken to be zero. Often, some nonzero elements of Y^"^^ and 
yCo!) dQ not occur in any complex invariants and hence do not occur on the right-hand side 
of ( p.29| ). These elements may be rephased freely. The Z matrix as defined in eq. ( p.30| ) 



has columns of zeroes corresponding to each of these elements. For simplicity, one may thus 
define a reduced vector Wr whose dimension is equal to the number of elements of Y^^\ 
f = u,d, which occur in complex invariants, and a corresponding reduced matrix Zr defined 
by 

^ = ZrWr (3.32) 

Clearly, by removing columns of zeroes we do not reduce the rank, so 

rank{Zr) = rank{Z) (3.33) 

Then the number of independent arguments among the complex invariants is given by the 
rank of Z: 

Nia = rank{Z) (3.34) 



The eqs. (|3.30| ) and (|3.32|) subsume the complex invariants of all orders (i.e. for the physical 



case Ng = 3, orders 4 and 6). One can also apply this method to determine the linearly 
independent invariants of a given order, 2n. For this purpose, we define a vector of complex 
invariants of this order: 

6n = (arg(X2n;i), ..., arg(Xjv_,2„))^ (3-35) 

We then have 

6n = Z2nW (3.36) 

As before, one can define a reduced vector Wr^2n composed of the elements of Y^^^ and Y^*^^ 
which occur in the complex invariants of order 2n, and a corresponding reduced matrix Zr,2n 
satisfying 

^2n = Zr^2nWr,2n (3.37) 

Then the number of linearly independent arguments among the complex invariants of order 
2n is given by rank{Z2n) = 'rO''nk{Z^^2n)- 

4 Further Results for General Nq 

In this section we give a number of new results on complex invariants and unremovable phases 
in Yukawa matrices for an arbitrary number of fermion generations, Nq- A consideration of 
general Nq yields further insight into the physical case Nq = 3. However, the results in this 
section will not be used directly in our analysis of realistic models, so that it may be skipped 
by the reader who is mainly interested in phenomenological applications. 

It is interesting to consider the case of Yukawa matrices with all arbitrary nonzero ele- 
ments, as we did before, in the context of our first theorem on the number of unremovable 
phases. Just as this case clearly yields the maximal number of unremovable phases, so also it 
yields the maximal number of complex invariants. We have proved above that in Corollary 
2 that in this case the complex quartic invariants form a complete set, in the sense that the 
arguments of the order 2n invariants for 2n > 6 can be expressed in terms of those of the 
quartic invariants. 
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We now calculate the maximal number of (independent) complex quartic invariants. This 
is done, as explained above, by performing the calculation in the case where the Y^-^^ have 
all nonzero arbitrary (complex) elements, since this case clearly yields the maximal number 
of such invariants. Our first result is that the maximal number of complex quartic P^^' 
invariants of each type f {= u or d) is 

{Np,j)^a. = (^^)' (4.1) 

Proof: 

For a given f = u oi d, the quartic P^^^ invariants are obtained by choosing all rectangles 
(some of which will be square). To pick a rectangle, we pick any two different columns of 
Y^'f^; this can be done in (2*^) ways. (Note that if we were to choose the two columns to 
coincide, we would get a doubly occupied line segment, corresponding to a real invariant.) 
These columns are equivalent (whence the division by 2!) because of the fact, noted above, 
that since each complex invariant and its complex conjugate yield the same phase constraint, 
we can count rectangles of the form || =|| and || =|| as equivalent. Within one of the 
columns we pick any two different elements; this can be done in any of (2 '') ways. Therefore, 
the number of (independent) complex quartic P'--^^ invariants of each type / for this case. 



which is the maximal number of such invariants, is as given in (|4.1| ). □ 

Since the maximal number of quartic P^-^-* complex invariants of each type f = u or d, 
Np^j^rnaxi IS the Same for f = u and f = d, the total maximal number of quartic P complex 
invariants is just double the right-hand side of ( |2.20| ), i.e. 

(A.,.)_ . »^i)]! (4,2) 

For the maximal number of complex quartic Q invariants, we find 

(A,,J_ . M^ (4.3) 

Proof: 

As above, the maximal value is calculated by considering the case where Y^'^^ and Y^'^'' have 

all nonzero arbitrary complex elements. First, we pick any column in, say, y("); this can be 
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done in A''^ ways. Within this column we choose two different elements, which can be done 
in (2^^) ways. This determines one side of the rectangle. To pick the other side, we choose 
any column of y^'^); this can be done in Nq ways. The two elements within this column are, 
of course, already determined by those in the column of F^"-*. As before, the cases || and 
It give complex conjugate invariants and hence the same phase constraint, so that one does 
not have to treat them together, which is equivalent to ignoring an overall reversal in the 
directions of all arrows. Hence, the number of complex quartic Q invariants is Ng(2'^)^g, 
i.e. the result given in ( [4.1|) . □ 

Hence, finally, the maximal value of the total number of independent complex quartic 
invariants is 

y-'-^ inv ,4 ) max \-^^P4)max ' y-'-^QiJmax 

N^{Ng-1){2Ng-1) 



2 
For the specific cases No =1, 2, and 3, we thus have 



(4.4) 



{{NpJ^ax, {NQ.Uax, (iVmM)™a:.) = (0,0,0), (2,4,6), (18,27,45) (4.5) 

Continuing with the analysis of this case in which all of the elements of Y'^^'> , f = u,d 
are nonzero arbitrary complex quantities, we immediately see that for this case (except for 
Ng = 1 for which {Ninv)max = {Np)max = 0) the resultant maximal number of unremovable 
phases, or, by our theorem on invariants, equivalently, the maximal number of indepen- 
dent arguments among the complex invariants, {Nia)max = {Np)maxi is always less than 
{Ninv,4)maxi SO that the arguments of the various complex quartic invariants are not all in- 
dependent. This is proved easily by recalling our earlier result (|2.2CI| ), which, with ( |2.17| ), 
yields 

{N,a)max = {Np^ax = (Ng - l)i2NG ~ 1) (4.6) 



Comparison with our result for {Ninv,4)max in eq. ( |4.4| ) shows that 

{N,nvA)max > {Nia)max f OT Ng > 2 (4.7) 
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We further recall from Corollary 2 to our theorem on invariants that in this case where all 
elements of the Y^^^ are nonzero arbitrary complex numbers, the arguments of the 6'th and 
higher order complex invariants can be expressed in terms of those of the quartic invariants. 
Although in this case the 6'th order invariants do not yield any further phase constraints, 
it is of interest to calculate the maximal number of these invariants since it shows how this 
number depends on Nq- 

For the number of independent 6'th order complex P invariants in this case of Y^-l"^ with 
all nonzero arbitrary (complex) entries, which gives the maximal number of complex Pg 
invariants, we get 

Proof: 

We calculate this as follows. First, we pick a column, in F*^-^-' (f = uoid), which can be done 
in Ng ways. Within this column, we pick two different elements, which can be done in (2 '^) 
ways. Next, we pick another column in this matrix, which can be done in (A^^ ~ 1) ways. In 
this column, we pick the element which is at the same level (i.e. has the same value of the 
first index) as the head of the arrow in the first column, or the element which is at the same 
level as the tail of this arrow; there is a factor of 2 corresponding to these two choices. We 
then pick a second element in this column, which is different from the first and is also not 
at the same level of the other end of the arrow in the first column (since if it were, we could 
immediately form a complex quartic invariant, i.e. we would not be constructing a genuine 
6'th order complex invariant). There are {Ng — 2) ways of choosing this second element. 
Finally, we choose the third column, which is different from the first two; this can be done 
in {Ng — 2) ways. The vertical levels of the elements in this column are now automatically 
determined from the previous choices. Since all of the columns are equivalent, we divide by 
3!. Putting these factors together gives, for the number of complex 6'th order P invariants, 
the result Ng{2''){Ng - 1) ■ 2 ■ {Ng - 2){Ng - 2)/3!, which yields (^. □ 
Including the f = u and f = d sectors, one obtains for the maximal value of the total number 
of Pq complex invariants just doubles the right-hand side of ( [4.8| ). 

For the maximal number of independent 6'th order complex Q'^^f^ ^ invariants, we calcu- 
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late 

(^,.,„,).... ^ ^°'^° - f ^° - ^' (4.9) 

Proo/; 

We consider the ¥^■1'^ matrix first. As before, we pick a column, and two different elements 
within this column, which can be done in Ng{2'^) ways. We then pick a different column 
within this matrix; this can be done in {Nq — 1) ways. In this column we pick the element 
which is at the same level (i.e. has the same value of the first index) as the head of the arrow 
in the first column, or the element which is at the same level as the tail of this arrow; there 
is a factor of 2 corresponding to these two choices. We then pick a second element in this 
column, which is different from the first and is also not at the same level of the other end of 
the arrow in the first column (since if it were, we could immediately form a complex quartic 
invariant, i.e. we would not be constructing a genuine 6'th order complex invariant). There 
are {Nq — 2) ways of choosing this second element. These two columns are equivalent, so 
we divide by a factor of 2!. Finally, we pick a column in the other matrix, Y^^ \ which can 
be done in Ng ways. The elements in this column are now completely determined by the 
previous choices. This yields the resuh A/'G(^'')(iVG - I) ■ 2 ■ {No - 2) ■ (2!)-^ ■ Nq, i.e. the 
number in ( [4.9| ). □ 

Including both the ///' = uud and ddu terms just doubles the number in ( [4.9|) . Hence, 
for the total number of independent complex 6'th order invariants in this case, we finally get 

{Ninvfi)max = 2[{Npgjjmax + i^QeJff'jmax] 

2N^{NG-lf{NG-2){2NG-l) 



(4.10) 



3 

Of course, in the case where the Y^^^ have all nonzero arbitrary complex elements, we have 
already shown that the quartic complex invariants form a complete set, so that these 6'th 
order invariants do not yield any new phase constraints. One use of the formulas ( [4.8| )-( ^rT0D 



is that they show that there are no complex 6'th order invariants for A^^ = 2. We will show 
below that all higher order complex invariants also vanish for Ng = 2. 
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The maximal number of independent complex invariants of order > 8 may be calculated 
in a similar manner. Here, as explained before, we consider only irreducible 8'th order 
invariants, i.e. those which cannot be factored into products of lower-order invariants. Thus, 
for example, 8'th order invariants formed as the products of two different quartic invariants 
are not counted. Since the proofs become rather long, we will simply give the results. For 
the maximal number of (irreducible) 8'th order P^-^-* invariants of each type f = u oi d we 
find 

[lyPsjJmax = ^j l4.iij 

For the maximal number of complex Qg ' invariants (which is, of course, equal to that of 
the Qg ' invariants), we calculate 

(N ) -(N ) _ NUNG-mNG-2){NG-3y 

(2 2) 

Finally, for the maximal number of complex Qg ' invariants, we get 

(^ ^ m{NG-mNG-2){NG-3) 

We now observe that the maximal numbers of each of these various types of 8'th order 
invariants, {Npgj)max in (|4.11|) , {N^(3,i))max = {N^(i,3))max in (|4.12| ), and {N^(2,2))max in 



^8 



( [4.13| ), all vanish unless Ng > 4. This explicit calculation thus provides an illustration of the 
general result in part (d) of our theorem on invariants. Recall that that result applied more 
generally to all complex invariants, including reducible ones, and stated that for the physical 
case Ng = 3 the arguments of any complex invariant of order > 8 could be expressed in 
terms of the arguments of complex invariants of order 4 and 6. Our general calculation above 
shows that, for the physical case A''^ = 3, there are no irreducible complex invariants of order 
8, i.e. all 8'th order complex invariants are products of lower-order invariants, and hence 
obviously their arguments can be expressed in terms of those of these lower-order complex 
invariants. 

One might be curious how part (d) would read if A''^ were equal to 2. The answer is: 
Theorem: 
For Ng = 2, the only (irreducible) complex invariants are of 4'th order. 
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Proof: 

To prove this, we calculate the maximal number of complex invariants as before, by using 

the case where the Y^f"* have all nonzero arbitrary complex elements. Consider first a P^ 



{/) 

2n 



complex invariant of order 2n > 6. To form this, one must pick n > 3 different columns from 
Y^f\ But this is impossible, since Y^^^ has only A''^ = 2 columns. Consider next a 6'th order 
complex Q-type invariant, say Q^"^'^^^ . To form this, one must pick two columns from Y^'^^ 
with the property that the ends of the arrows do not match each other (if they did, one would 
form a P4 invariant). But this is impossible since F'-"-* is just a 2 x 2 matrix. Similar reasoning 
shows that it is also impossible to form any higher-order irreducible complex invariant. □ 
Of course, there may be reducible complex invariants of higher order, e.g., products of two 
different quartic invariants. 

For this case Nq = 2 it is easy to construct all such invariants. To obtain the maximal 
set, we assume that the Y^f'^ have all nonzero arbitrary (complex) elements; obviously, if 
some elements are zero, the resultant set of complex invariants would be commensurately 
reduced. According to our general result (|4.1|), {Np_^j)max = 1, and the explicit complex 
invariants are 

-'11,22 — -f 11 -'^22 -'21 -^12 l^'-L^J 

for f = u,d. From our general result (|43| ) it follows that there are {NQ^)max = 4 independent 
complex quartic Q invariants. These are 

Qii,22 = Ylr^Y,fYir^*Y}^* (4.15) 

gi2,2i = Yl-^Y^?Y,^tY}?* (4.16) 

Qii,2i = Ylr^Yit^Yir^*Ylt^* (4.17) 

and 

V12,22 — -'12 -'22 -'22 -^12 l^--'-OJ 

From eq. (|2.2CI|) , we know that there are A'p = 3 unremovable phases for this case, and 
correspondingly, among the Ninv = 6 complex (quartic) invariants there are Nia = Np = 3 
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invariants with linearly independent arguments. These may, for example, be taken to be 

-^11,22) -^11,22) aild (5ll,22- 

Returning to the case of general Nq, as one makes some of the elements of the Y^"^^ 
and Y^'^' zero, the number of unremovable phases and the number of complex invariants 
are decreased. The amount by which the number of unremovable phases is decreased is 
not solely a function of the number of zeroes, but in addition depends on their locations. 
For almost all of the models which we have studied, one finds a simple empirical relation 
Np = {Np)max — Nz = {2Ng — 1){Ng — 1) — Nz, where we recall that Nz is the total number 
of zero elements in Y^"^^ and Y'^'^\ This is not a general result; in section 5 we exhibit a 
toy model (no. 8) which, considered as an Nq = 3 model, is an exception. However, that 
model is degenerate, in the sense that it involves a decoupled third generation in both up 
and down quark sectors. If one considers only the mutually coupled sector of the model, 
then it reduces to a Nq = 2 model, and does obey the Nq = 2 version of the above empirical 
relation. 

As a final item in this section, we contrast the unremovable and hence physically mean- 
ingful phases in the quark Yukawa (or mass) matrices, with the unremovable phase (s) in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix, V, defined by the charged weak 
current 

J^ = UjL,ml^VjkdkL,m (4.19) 

where UjL,m and dkL,m are the left-handed components of the mass eigenstates of the Q = 2/3 
and Q = —1/3 quark fields (indicated by the subscript m) defined after the spontaneous 
symmetry breaking of SU{2) x f/(l), with UiL,m = UL,ni, U2L,m = Cm, M3L,m = tL,m, and 
similarly for dj^m- If the mass matrices are diagonalized by the biunitary transformations 
Uf,LY^'^'^Uf^j^ = Y^lJg for f = u,d, then V = Uu,lUIl- The rephasing properties of V are quite 
different from those of the Y^'^', as is obvious from the fact that the charged weak current is 
chirality-preserving, connecting left-handed with left-handed chiral components of individual 
physical mass eigenstates of quark fields, whereas the Yukawa terms connect left-handed 
with right-handed chiral components. Furthermore, the Yukawa matrices, being defined via 
SU{3) X SU{2) X t/(l)-invariant Yukawa couplings, involve the actual SU{3) x SU{2) x U{1) 
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fields, i.e. the Nq left-handed SU(2) doublets Qjl and the 2Nc right-handed SU(2) singlets 
UjR and djn (not mass eigenstates) . The rephasing equations for the Yukawa matrices were 
given above in ( p.5| )-( p77D . In contrast, when one determines the unremovable phases in V, 
since it is defined in terms of the mass eigenstates themselves, the rephasing equations which 
one uses are 



j,Lm — ^ ' "'j,Lm 



^j,Lm € U 



d,,Lm = e^^^c?;.^^ (4.20) 

for j = 1, ..., Nq. The number of these rephasing equations is 2Nc, as opposed to the 3A^g 
rephasing equations ( |2.5|) -( pr7D for the Yukawa matrices. The resultant CKM matrix is 

V^^ = e'^-'^+^'^^V.k (4.21) 

Recall the well-known argument that since (i) V is unitary, and hence specified by A''^ real 
parameters, (ii) there are 2Ng rephasings as noted; and (iii) the overall rephasing U(l) 
defined by ( [4.20|) with 9j = (pj for all j = 1, ■■■,Ng leaves V invariant, it follows that there 
are Nq — {2Ng — 1) = {Nq — 1)^ real parameters specifying V , of which Ng{Ng + l)/2 are 
rotation angles and the remainder are 

N,, ^ '^° - ^f ° - ^> (4.22) 

unremovable, physically meaningful, phases. An immediate contrast is that in general, Npy 
is only a function of Ng, whereas the number of unremovable phases in the Yukawa matrices 
do not just depend on Nq but rather on the detailed structure of these matrices. An 
important inequality is that the number of unremovable phases in the quark Yukawa matrices 
is at least as large as the number of unremovable phases in the CKM matrix V: 

Np > Npy (4.23) 

The proof of this is obvious since the phases in V originated, as V itself did, from the diag- 
onalization of the Yukawa (equivalently, mass) matrices as defined in ( |2.1|) . Our discussion 



27 



of specific realistic models will illustrate this inequality; in particular, in the physically in- 
teresting case Ng = 3, the CKM matrix has Npy = 1 unremovable (CP- violating) phase, 
whereas the realistic models typically have at least Np = 2 unremovable phases in the Yukawa 
matrices. From our inequality ([4.23|) , there follows the theorem that phases in quark mass 
matrices do not necessarily violate CP. The issue of CP violation involving Yukawa couplings 
also depends on how many Higgs are present in the theory and how they couple to the quarks. 
In order to make our analysis as general as possible, we have defined the mass and Yukawa 
matrices as in eq. ( |2.1|) without explicitly specifying the precise Higgs content of the the- 
ory; however, most recent studies of actual models (including our own) assume the minimal 
super symmetric standard model, for which the Higgs were listed in ( ^.31 ) and ( ^.41 ). For a 
typical case where Np > Npy, i.e. Np > 1 for the Nq = 3 case of physical interest, it follows 
that some phases in the quark mass matrices may contribute to CP-conserving quantities. 
In the mathematical case No = 1, eqs. (|2.2CI| ) and (^4.22|) imply that {Np)max = Npy = 0. 



For the case Nq = 2, the CKM matrix still has no phase, but ( ^.2(j| ) shows that Np may be as 
large as 3. Indeed, in the specific Nq = 2 case worked out completely in eqs. ( [4.14| )-( ^rT^ ), 
Np is equal to 3. Thus, in this A''^ = 2 case, the phases in the Yukawa (mass) matrices do 
not yield any phase in the CKM matrix. 

5 Applications to Specific Models 

Although our theorems are quite general, it is useful to see how they apply to various specific 
models. We shall do this in the present section. Over the years, a large number of mod- 
els have been proposed for the origin of quark masses and mixing. In some of the earliest 
efforts, the forms of the models were hypothesized to apply near to the electroweak mass 
scale. In general, various studies can be classified according to whether they assume a the- 
oretical framework of perturbative or nonperturbative electroweak symmetry breaking. We 
shall concentrate here on the class of models in which the observed electroweak symmetry 
breaking is perturbative. In this class, an appealing framework is provided by supersymmet- 
ric extensions of the standard model, which stabilize the Higgs sector and hence stabilize 
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the hierarchy between the electroweak scale and the Planck scale. Among these supersym- 
metric extensions, the minimal supersymmetric standard model (MSSM) has the particular 
feature that it yields gauge coupling unification, thereby providing a self-consistent basis for 
a supersymmetric grand unified theory (SUSY GUT) |Q. In turn, for studies of fermion 
masses and mixing, a (supersymmetric) grand unified theory has the advantages that (i) 
one can naturally relate quark and lepton masses; (ii) with the grand unified group SO(IO) 
one can naturally obtain (complex) symmetric Yukawa matrices at the GUT scale, which 
helps to minimize parameters. Thus, a number of studies have been carried out of particu- 
lar (complex) symmetric forms for Yukawa matrices at the (supersymmetric) grand unified 
scale, in which the renormalization group equations of the MSSM are used to evolve these 
forms down to the electroweak scale, where they are diagonalized to yield the known quark 
masses and the CKM quark mixing matrix V. In these studies, one assumes some sym- 
metries (which are usually discrete in current work) to prevent various Yukawa couplings 
and thereby render various entries in the Yukawa matrices zero. The purpose of this is, of 
course, to minimize the number of parameters and hence increase the apparent predictive- 
ness. (Parenthetically, we recall that when one evolves these forms for the Yukawa matrices 
down to the electroweak level, they do not, in general, remain precisely symmetric, and some 
of the zero entries can become nonzero.) These studies are phenomenological in the sense 
that they do not explain the origin of these discrete symmetries. The most likely origin of 
the discrete symmetries is probably an underlying string theory (in which these are actually 
local symmetries). However, ironically, string theories do not generically yield either simple 
gauge groups (the essence of grand unification) or symmetric (complex or real) mass matri- 
ces, at least at the present limited level of understanding of their physical properties P . At 
a phenomenological level it is reasonable to seek forms for the F'--^^ with as many zero entries 
as possible. However, if and when one understands fermion masses and mixing in terms of 
a truly fundamental theory, with perhaps no free parameters, then since everything is (at 
least in principle) calculable, it makes no difference whether the elements of the Y^^^ are 
exactly zero or just small. The entries which are modelled as being exactly zero might well 
be nonzero, calculable quantities suppressed by sufficiently high powers of (Mr/Mp), where 
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Mr is the reference scale at which one analyzes the effective Yukawa interaction. Thus, the 
real goal is not necessarily to maximize the number of zero entries in the Yukawa matrices 
of a viable model as one necessarily tries to in a phenomenological approach; rather, it is 
to gain an understanding of the (presumably unique) fundamental theory, together with a 
calculational ability which would enable one to calculate these matrices from first principles. 
Given the widely acknowledged [Q limitations on one's understanding of the phenomeno- 
logical consequences of strings and, at the level of this understanding, the apparent lack of 
uniqueness of string theory, as evident, e.g., in the 4D free-fermion constructions, it does not 
appear that one is very close to achieving this goal at present. Finally, for completeness, we 
must note that there are a number of concerns which must be faced in actual model-building 
with super symmetric theories, such as, for example, the "/i problem" and the problem of 
splitting light Higgs doublets from other components of grand unified Higgs representations 
which must be heavy, and so forth; these are reviewed, e.g., in |^. 

With these brief theoretical comments, we proceed to illustrate our general theorems with 
specific models. We emphasize that our purpose here is not to propose new viable models of 
quark Yukawa matrices with the attendant renormalization group calculations and compari- 
son with experiment; rather, our purpose is to consider (generalizations of) forms which are 
known to be experimentally viable and, for these, to determine the unremovable phases and 
which elements can be made real by rephasings. Since our results do not depend on whether 
or not the Yukawa matrices are symmetric at some mass scale, we shall consider the general 
case of non-symmetric Yukawa matrices.^ As a basis for this, we will use a valuable recent 
study of SUSY GUT Yukawa matrices with the restriction that |l^fc | = 1^^^ | [Q- This work 
presented an exhaustive list of five viable models, with this restriction, which we label Mj, 

j = 1 — 5. For the reasons mentioned above, we shall actually consider the generalizations 

.(f). \ ( f) \ 

of these models in which 11^^! is not necessarily equal to ll^j |. Again, we emphasize that 

this does not affect the counting of the number of unremovable phases or the determination 



^From a SUSY GUT model-builder's point of view, this might seem perverse, since symmetric Yukawa 
matrices have the appeal of reducing the number of parameters and, moreover, can be naturally obtained 
in grand unified theories. But this appeal is offset by one's awareness of the results of string-based studies, 
mentioned above, which suggest that nature may not, in fact, be described by either a grand unified group 
or by symmetric Yukawa matrices. 
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of which elements can be made real by rephasing. We label our five generalizations as M', 
j = 1 — 5 and discuss them in the next five subsections. Since the respective special cases are 
experimentally viable, it follows, a fortiori, that the generalizations are also viable, i.e., for 
appropriate choices of the parameters, one finds predictions in agreement with experimental 
constraints on the CKM matrix for experimentally acceptable input values of the top quark 
mass. Our set of models is not intended to be exhaustive, since a comprehensive study 
of viable non-symmetric Yukawa matrices has not been performed (for the understandable 
reason that at a phenomenological level such models have too many parameters to make 
useful predictions). Model 6 is related to a model with non-symmetric Yukawa matrices re- 
cently sketched in a string context [^. Models 7 and 8 are toy examples to illustrate certain 
theoretical points. 

5.1 Model 1 

The first illustrative model is defined by the Yukawa matrices 

Ai2 
r(") = I A21 A22 I (5.1.1) 

A33 

/ fii2 \ 

F(^) = B21 B22 B23 (5.1.2) 

\ B32 -B33 / 
where here, and below, each element in the Yukawa matrices is, in general, complex, since 

no symmetry requires it to be real. This model has Nz^u = 5, A^^^^ = 3, A^^ = 8, and 

iVeq = 10. In the physical context, the forms in eqs. ( |5.1.1| ) and ( p.l.2| ) are assumed to hold 



at a given mass scale, such as the SUSY GUT mass scale Mqut = (1 " 2) x 10^^ GeV where 
the gauge couplings coincide in the simple MSSM calculation, or a higher mass scale near to 
the (reduced) Planck mass Mp = 2.4 x 10^^ GeV, which would be expected in supergravity 
and string theories, in which one would extend the MSSM calculation to take account of 
string threshold corrections and additional fields which could be present. Special cases of 
this model have been studied in Refs. and [§]. Besides the assumption made in both of 
these works, that \YJP\ = \Y^p\, f = u,d, Ref. also took l^asl = l^ag]. 
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We will give a detailed, pedagogical discussion of the application of our methods for this 
model. We first calculate the T matrix. The reader may recall here the definition in eq. 
( p.l3| ) and the fact that the ordering of the columns of T is determined by the definition of 
the vector v in (|2.10|) while the ordering of the equations in (p.9|) , and hence the ordering 
of the elements of the vector w and the rows of T, is defined by the sequence of nonzero 
elements Y^P in J^T^ } with jk = {11, 12, 13, 21, 22, 23, 31, 32, 33} for / = u and then f = d. 



Ti 



(5.1.3) 



Explicitly, here the rows correspond to the rephasing equations for Y12 , Y21 , Y22 , ^33 , 
Yif, ^2?\ ySK Y23\ Yaf, and Y^^l We obtain the iV,, x 3Ng = 10-row by 9-column matrix 

/100010000\ 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 
10 1 
10 10 

Vooioooooiy 

where the subscript refers to the model. Further, we calculate the rank of this matrix, 
rank{Ti) = 8. We can then apply our theorem in ( p. 17] ) to conclude that this model has 
Np = Neq — rank{T) = 2 unremovable, and hence physically meaningful, phases in the 
Yukawa matrices Y^-^^ (and thus of course also in the mass matrices). According to our 
theorem on invariants, there are therefore Nia = Np = 2 complex invariants with independent 
arguments. For this model we find that A^m„,4 = N^a, i.e. the number of quartic complex 
invariants is equal to the number of complex invariants with independent arguments, so 
that the two complex quartic invariants fully account for the two unremovable phases and 
corresponding phase constraints. Explicitly, we find the (nonzero, independent) complex 
quartic invariants 



P. 



id) 



22,33 



{d)^{d)^{f)*^{f)* 

-'^22 -'33 -'^32 -'^23 



-022-033-032-^23 



(5.1.4) 
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and 



V12,22 — -'12 -'22 -'^22 ■'l 



12 

= AuB22a;^bi^ (5.1.5) 

The graphical representations of these complex invariants were given above in Fig. 1(a) and 
Fig. 1(b), respectively. Associated with these are their two rephasing- invariant arguments 
(phases), arg(P2233) ^iid arg((5i2,22)-0 These complex invariants and associated invariant 
arguments imply constraints on which elements of Y^"^^ and Y^'^' can be made real. Since, in 
general, arg(P22,33) ¥" 0) ^) it follows that (i) at least one of the Np = 2 unremovable phases 
must reside among the set 

V ={^22, 523,532,^33} (5.1.6) 

-^22,33 

and (ii) the 2x2 submatrix in Y^'^^ formed by the set ( 5.1.6| ), 



\ -D32 -D33 / 

and hence also Y^'^'^ itself, cannot be made real or hermitian.[^ Second, since, in general, 
arg((5i2,22) 7^ 0, 71, it follows that the elements in the set 

'S'qi2,22 = {^12,^22,512,-822} (5.1.8) 

(iii) cannot in general be made simultaneously real and thus, (iv) if one chooses Y^^^ real, 
then it is not possible to make B12 and ^22 both real. These constitute the complete set 

^In the model of Ref. |^ the phase Oq^^^ 22 — arg(Qi2,22) did not significantly help the fit to experimental 
data and hence was taken to be zero to minimize the number of parameters. 

^° As noted in Ref. H, our phase results disagree with those in Ref. m for the five models which they 
consider. In particular, we disagree with (i) their statement (footnote on p. 21 of Ref. [^) that (complex) 
3x3 symmetric mass or Yukawa matrices can be transformed to hermitian form by rephasings of the quark 
fields; (ii) the counting of unremovable phases in the Yukawa matrices in that work; and (iii) the claims 
made concerning which elements of these matrices can be made real by rephasings of fermion fields. For the 
models Mj, j = 1 — 5 considered here, and hence also for the corresponding models Mj of Ref. 0, we find 
the values of Np are respectively 2, 3, 2, 2, 2, not 1, 2, 1, 1, 1 as claimed there. In Ref. the rephased forms 
of the Yukawa matrices (their Table 1) list F*^''^ as being hermitian and F*^"^ as being real for all five models. 



We find that Y^'^^ cannot in general be made hermitian (nor can the submatrix ( 5.1.7 ) be made real or 
hermitian) in models 1, 2, and 3. We also find that y*-"-* cannot be made real (or hermitian) in models 4 and 
5. We have communicated our results to P. Ramond and G. G. Ross and thank them for discussions. Our 
differences do not invalidate the important conclusions of Ref. Q] that these five forms are experimentally 
viable, since the inclusion of more phases can only improve the fit to experiment. 
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of rephasing constraints on the Y^-l'\ f = u,d. These constraints allow both phases to be 

put in Y^'^^ and both to be put in the set Sp(d) . If one chooses to make B22 and B12 real, 

22,33 

then one cannot make Y^^^ real, and must assign one phase to A12 or A22 and the second to 
B23, -B32 or B33. The elements A21, A33, and B21 do not occur in any complex invariants; 
consequently, their phases are unconstrained and may be rephased to arbitrary values by the 
rephasing transformations ( |2.5| )-( P77| ). Our results show that it is not, in general, true that 
a (complex) symmetric Yukawa or mass matrix can be transformed to a hermitian matrix 
by the rephasings of fermion fields. 

As an explicit illustration, after rephasings, one could obtain 



yW 



/ 

l^2l| 





IA2I 
1^22! 








I ^33 1 



(5.1.9) 



and 



yW 



( 

I-B21I 

V 



\B12\e 
\B22\e 



-i arg(Qi2,22) 
i arg(P2<2'33) 



IB 



32 1 



\ 

I-B23I 
l-^ssl / 



(5.1.10) 



where for any real quantity, \Y!>1'\, one could have obtained 



jk h 



Y^l \ due to the freedom of 



Af) 



If) 



W^ 



using rijl = 1 rather than rijl = in ( p.9| ). Since arg(Y2i ) = arg(i?2i) is unconstrained, one 



could also have obtained 1^21 = |-B2i|e* '^''s ('212,22) [y^ gq_ ( |5,1,10|) ^ as discussed in connection 



with eq. (|2.18|) , although of course this choice has no effect on the physics. An example of 
a form which could not, in general, be obtained from any rephasings of fermion fields is the 
hermitian form 

yW ^ 




I R l^«5 

|-Di2|e 
I-B22I 



(5.1.11) 





|-D23|e ^ 

l-B32|e~*'^ I533I 

for any values of 6 and (j). As is evident from our invariants, this statement does not depend 
on whether or not the initial y-"-* and Y^'^^ are (complex) symmetric. 

For this model, we shall also illustrate explicitly the operations which enter into the 
proof of our theorem ( ^.17| ). In order to get the invertible matrix T, we select and delete 
Npjj — r^ = 10 — 8 = 2 rows from the T matrix for this model, given in eq. (|5.1.3| ). 



<eq 



This corresponds to deleting two elements of the vector w defined in ( 2.11 ) to get w, i.e. 
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not rephasing the corresponding Yj^' . Given that we have constructed the full set of Nia 
invariants with independent phases, we know the full set of allowed choices for elements to 
be rephased and correspondingly the choices for the Np = 2 elements which will not be 
rephased. In this model, the rephasing equations for the elements comprising the complex 
invariants ^22,33 ^^^ ^12,22 correspond to the rows (1,3,5,7) and (7,8,9,10) of the matrix T, 
respectively. Hence we must delete one row from the first set and one row from the second 
set; i.e. we have the following 15 choices of row pairs to delete: (rl,r2), with rl chosen 
from (1,3,5,7) and r2 chosen from (7,8,9,10), excluding, of course, the choice rl = r2 = 7. 
For our specific illustration, we delete rows 5 and 7 and thus do not rephase Y12 or 1^22 ■ 
Note that if one were to try to remove two rows in a manner which would violate the phase 
constraints, one's error would immediately be signalled by a reduction in the rank of the 
resultant matrix. For example, if one were to try to remove rows 1 and 2, one would obtain a 
matrix with rank 7. Recall that this choice is forbidden since it would entail the implication 
that one could rephase all other elements of Y^'^^ and Y^'^\ which is false, since it would mean 
that one could in general rephase the argument of Qi2,22 to zero, contrary to fact. Note that, 
in principle, one may determine which rows may be deleted simply by testing the rank of the 
resultant matrix, without first constructing the complete set of complex invariants. However, 
obviously the more efficient procedure is to construct this set of invariants first, since once 
this is done, it is clear which sets of Np rows can or cannot be deleted consistent with the 
phase constraints. Continuing with our explicit calculational example, having deleted the 
rows 5 and 7 from T, and thereby obtained an r^ x 3Ng = 8-row by 9-column matrix, 
we next delete 3Ng — r^ = 1 column from this matrix in order finally to get the square, 
Tj, X Ty = 8 X 8 matrix T. This corresponds to deleting one element of the vector v defined 
in eq. ( p. 10 ) to get v. We find that any of the nine columns of T may be deleted in this 
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step; for definiteness, we pick the ninth. This yields the reduced matrix T for this modeh 



fi 



/10001000\ 

10 10 

10 10 

10 1 

10 1 

10 

10 1 

yooiooooo; 

(with determinant —1). By construction, this is invertible, and we find 



(5.1.12) 
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(5.1.13) 



As an illustration of our theorem on invariants, we discuss the 6'th order complex invari- 
ants for this model. Since we have already constructed as many independent complex invari- 
ants with independent arguments as there are unremovable phases, using quartic complex 
invariants, our theorem on invariants and its specific Corollary 1 imply that the argument of 
any 6'th, as well as higher order, complex invariant, are expressible in terms of the arguments 
of the quartic invariants. By explicit computation, we find that there is one complex 6'th 
order invariant (i.e. Ninvfi = 1): 



V32,23,12 



-832-823^12-833-812^22 



(5.1.14) 



Thus the total number of independent complex invariants for this model is Ninv = Ninv,A + 
Ninvfi = 3. The graphical representation of this complex invariant is shown in Fig. 3(a) 
(where here and elsewhere, lines which cross each other at a point where there is no or 
®, as in the Y22 position, do not mean that the invariant contains this element). Note that 
this 6'th order invariant satisfies the condition (i) of part (d) of the above theorem. It is 
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thus easy to find the exphcit reduction; we obtain the factorization 

arg(532523Ai2i?3*35i2^;2) = ^Tg{B^2B2zA^2{B22Bl^)Bl^Bl^Al^) 

= aTgiA,2B22A;^Bl^) + aigiBs2B23B;^B;.,) (5.1.15) 

whence 

arg((53?23,i2) = arg((5i2,22) - arg(P2^2?33) (5.1.16) 

Indeed, in this case, we find the factorization 

^32,23, 12 = 1-^22 1 ^12,22-^23,32 (5.1.17) 

Recall here that P2332 — -^223*3- Note also that according to our definition given above, 
Q32J3 12 is still an irreducible 6'th order complex invariant because it is not expressible only 
as a product of lower-order complex invariants; there is also the |-B22|~^ factor. Graphically, 
the insertion of the factor B22B22 in eq. (|5.1.15|) is equivalent to inserting a coincident arrow 



head and tail, or and ® in the 1^22 position. Having done this, one sees immediately that 
the 6'th order invariant factorizes in the manner specified by eq. ( |5.1.17| ). In addition to 



using our theorem on invariants to show that the argument of Q\,2 23 12 can be expressed in 
terms of those of quartic complex invariants, another way to see this is to use the Z matrix 
method discussed above. For this, we first define the vector ^ defined in eq. ( p. 311 ) for the 
present model: 

(Oi = (arg(P2?33), arg(Qi2,22), arg(QgSi2))^ (S-l-lS) 

where the subscript on (,^)i denotes the fact that this is the ^ vector for model 1. The 
ordering of the elements in the vector w was given in general after its definition, eq. ( |2.11| ), 



and explicitly for this model. Then Z is an Nin^ x N^.^ = 3-row by 10-column matrix. From 
the definition ( |3.30| ) we calculate it to be 
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-1 


-1 1 
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-1 


-1 1 
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-1 


-10 
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1 1 



(Z)i =10-10-1010 (5.1.19) 
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where again the subscript on {Z)i refers to model 1. Since the elements A21, A33, and B21 
do not occur in any complex invariants, they may be removed from the vector w to yield 
Wr, and correspondingly, one may remove the second, fourth, and sixth columns from Zi to 
form the reduced 4-row by 7-column matrix Z^ defined in ( p.32| ). 



(Z; 



r 1 



1-1-11 
1-1-110 
1-1-10 1 11 



(5.1.20) 



Since we have only removed columns of zeroes, this does not change the rank of the matrix. 
We calculate that rank{{Z)i) = rank{{Zr)i) = 2, in agreement with our theorem that 
Np = Nia and the relation A^^ = rank{Z) from eq. (|3.34|) . We find further that removing 
the last row from {Z)i or {Zr)i does not reduce the rank, which shows that the argument 
of Q32,23,i2 can be expressed in terms of those of ^22,33 ^^d (5i2,22- The phase properties 
and complex invariants of the Yukawa matrices for this and the other realistic models to be 
discussed are summarized in Table 1. Having given a detailed discussion of the details of our 
analysis for this model, we next proceed with the others, which we will discuss more briefly. 



5.2 Model 2 



This model is defined by 



y(") 



and Y^'^^ as in model 1, i.e. 



yW 



/ A12 

A21 A23 
\ A32 A33 



B12 

B21 B22 B23 

B32 -B33 



(5.2.1) 



(5.2.2) 
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model 


Ne, 


^T 


iVp = AT,, 


^^inv,A 


^4 


^4 


Ninv,& 


Pe or Qe 


1 


10 


8 


2 


2 


-'22,33 


<5l2,22 


1 


^{ddu) 
V32,23,12 


2 


11 


8 


3 


4 


-'22,33 


<5l2,32, <523,32, ^23,33 


2 


V32,23,12? V22,33,12 


3 


10 


8 


2 


2 


-'2233 


<5l3,32 


1 


^(ddu) 
V22,33,13 


4 


10 


8 


2 


2 


pH 

-^^22,33 


Ql2,22 


1 


V32,23,12 


5 


10 


8 


2 


2 


p{«) 
-* 22,33 


Ql3,22 


1 


^{uud) 
V32, 13,22 


6 


12 


8 


4 


7 


p{«) 
-'^22,33 

p{d) 
-'22,33 


Ql2,22, Ql2,32, ^22,32, 
Q22,33! Q23,33 


4 


V22,33,125 V33,12,22 

^ (ddu) ^ (ddu) 
V22,33,12) V33,12,22 



Table 1: Summary of unremovable phases and complex invariants for the six realistic models 
of Yukawa matrices. 



39 



This model has Nz^u 
matrix 



4, AT, 



z,d 



3, and hence A^, 



eg 



11. We calculate that T is the 11 x 9 



To 



/100010000\ 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 
10 1 
10 10 

Vooioooooiy 



(5.2.3) 



with rank 8. Hence, from eq. ( p . 1 7| ) , there are Np = N^g — rank{T) = 3 unremovable phases 
in the Y^'f\ f = u,d. By part (c) of our theorem on invariants, there are correspondingly 



Ni, 



N^ 



3 independent phases of complex invariants. We actually find that there are 



Ninv,4 = 4 independent complex quartic invariants in this model. (Recall from eq. ( |3.28| ) that 
in general there may be more independent complex invariants than independent arguments 
of invariants.) The four independent complex quartic invariants are P22 33 ^^ i^ model 1, 
together with 

Qi2,32 = AuBs2A;^BI^ (5.2.4) 



Q 



23,32 



and 



Q 



23,33 



^23-832^33-822 



^23-833^33-823 



(5.2.5) 



(5.2.6) 



Of these, the arguments of the first, third, and fourth invariants are related according to 



arg(P2?33) + arg(Q23,32) - arg(Q23,. 



33 J 







(5.2.7) 



so that, as required by our theorem, the number of independent arguments of invariants is 
Nia = 3. One can also determine the number Nia without using part (c) of our theorem and 



Corollary 1 on invariants by calculating the matrix Z defined in eq. ( |3.30|) and computing 
its rank. In the present case, it is sufficient, and simpler, to deal with the matrix Z4 defined 
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in eq. ( p.36| ) for quartic complex invariants. As defined in eq. (p.35|) for 2n = 4, the vector 



^4 of arguments of quartic invariants has dimension Ninv,A = 4 and is 

({4)1 = (arg(Pi2,33), arg(Qi2,32), arg(Q23,32), arg(Q23,33))^ (5.2.8) 



where the subscript 2 on (^4)2 refers to model 2. Using ( |3.36| ) and ( |2.11| ) with the specified 



ordering of elements in W4 (and rjjjj = as noted before), we then calculate the Ninv,4 x N^q 
= 4-row by 11-column Z4 matrix as 



iZ. 



4 2 



f ° 

















1 


-1 


-1 1 \ 


1 








-1 





-1 








1 








1 





-1 





-1 





1 


\o 





1 





-1 








-1 


1; 



(5.2.9) 



This matrix has rank 3, so that the three independent arguments among the four quar- 
tic complex invariants account for the full set of three unremovable phases in this model. 
By selecting and deleting rows of (.^4)2, one easily establishes which of these quartic in- 
variants have linearly dependent arguments; deleting row 1, 3, or 4 does not reduce the 
rank, while deleting row 2 yields a matrix with rank 2. This means that any one of the 

sets ((5l2,32, ^23,32, ^23,33), (A?33> Ql2,32, Q23,33), and (Pia,^) Ql2,32, ^23,32) COUStitutCS a 

complete set of Nia = 3 complex invariants with linearly independent arguments, while the 
arguments of the invariants in the set (-P22,33' ^23,32, <523,33) are linearly dependent, in agree- 
ment with the explicit calculation ( [5.2.7| ). We will take the three independent unremovable 
phases in this model to be arg(P22,33)5 arg((5i2,32), and arg((523,32)- Since the elements A21 
and B21 do not occur in any complex invariants, they may be rephased freely, and one may 
also form the vector Wr and corresponding matrix Zr^^ which is just Z4 with the second and 
seventh columns of zeroes removed. The constraints on rephasing from P22 33 ^^^ (i) ^^^ (ii) 
as in model 1, and (iii) none of the sets 

SQr,,s, = {^12, A32, P12, P32} (5.2.10) 

SQ23,32 = {^23,^33,-822,-832} (5.2.11) 

and 

'^Q23,33 = {^23,^33,523,533} (5.2.12) 
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can be made simultaneously real. In particular, if Y^'^^ is made real, then none of the sets 
{512,-832}, {-B325-B22}, and {523,-833} can be made simultaneously real. As in model 1, 
because of the constraint from P22 33' ^'''^^ cannot in general be made real or hermitian. An 
example of a rephasing of the Yukawa matrices allowed by these constraints is y*-"-* real, 

F(«)' = I IA21I IA23I I (5.2.13) 




and 

/ |5i2|e-*^''sWi2.32) \ 

1^21 1 |522|e-*"'"sW23,32) 1^231 



yW 



\ I-B32I |533|e*(^''S(^22,33)+arg(Q23,32)) I 



(5.2.14) 



The construction of T from T is performed in a manner similar to that discussed for model 
1. 

As in model 1, since we have already constructed as many independent complex invariants 
with independent arguments as there are unremovable phases, Nia = Np = 3, using quartic 
complex invariants, our theorem on invariants and its Corollary 1 imply that the arguments 
of any 6'th, as well as higher order, complex invariants, are expressible in terms of the 
arguments of the quartic invariants and hence yield no new phase constraints. We find that 
in the present model there are two different 6'th order complex invariants (so A^j„t),6 = 2): 

^32,23,12 = ^32^23-812^33^12-^22 (5.2.15) 

and 

^22,33, 12 = -S22 -833^12-^23 -^12^32 (5.2.16) 

The total number of independent complex invariants in this model is therefore Ninv = Ninv,4+ 
Ninvfi = 6. The explicit expressions for the arguments of these 6'th order complex invariants, 
in terms of quartic complex invariants, are 

arg((332!23,i2) = arg(g23,32) - arg(gi2,32) (5.2.17) 

and 

arg((3Siti2) = ^T^&iHtss) + arg(gi2,32) (5.2.18) 

42 



Again, a different way of sliowing tliat tlie 6'tli order complex invariants yield no new 
constraints is by analyzing the full Z matrix for the model. The method should be clear 
from our previous example, so we omit the details. 

Parenthetically, we consider the special case of this model when 1^22 ~ -^22 = 0. Here 
y *^") and Y'^'^^ have the same form. Although current interest is in models of Yukawa matrices 
which apply at the SUSY GUT level, it should be noted that the further specialization in 
which one takes IVj^! | = l^fci | is formally the same as one of the early works on quark mass 



matrices [|T0| . Of course in that work the form was hypothesized to apply near the electroweak 
level in a (non-supersymmetric) SU{2)l x 5'f/(2)/j x f/(l) electroweak gauge model, not a 
SUSY GUT. As is well known, whether one considers the symmetric form with B22 = to 
apply at the electroweak level or at the SUSY GUT level with MSSM evolution equations, 
it is currently disfavored by experiment, because of the large value which it predicts for Vet, 
given the lower bound on top quark mass indicated by both direct limits from CDF and DO 
rT| and by precision fits to LEP data [|r^. It is, however, of historical interest. In this case. 



Neq = 10, rank{T) = 8, whence Np = 2. Of the four quartic invariants for model 2, two 
vanish for this special case: P22 33 ~ ^^^ Q23,32 = 0. The other quartic complex invariants, 
<5i3,32 and ^23,33, remain, so that Ninv,i = 2. The phases of these are independent, so these 
two quartic invariants account for both of the two unremovable phases and resultant phase 
constraints. The single 6'th order complex invariant for model 2 vanishes for this special 
case, which therefore has no nonzero complex 6'th order invariants. 

5.3 Model 3 

The third model is given by 

An 
yH = I A22 I (5.3.1) 

^31 A33 



with Y^'^'> as in model 1: 



5i2 
yW = I S21 B22 B23 I (5.3.2) 

-B32 -B33 
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\Y^j \ were studied in Refs. JT^ and [^; in Ref. 



( f) \ 
Special cases of this model for which \Yjl \ 

|T^, the further assumptions IA22I = l^isl and I-B23I = 2I-B22I were made. The calculation of 

T should be clear from our previous examples, so we just mention that it is a 10 x 9 matrix 

and we find its rank to be 8. Hence, there are Np = 2 unremovable phases in the Y^^\ 

f = u,d. We find Ninv,i independent complex quartic invariants; these are ^22,33 ^^ i^ the 

previous models, and 

Ql3,32 = ^13-832^33-812 (5.3.3) 

These two complex quartic invariants clearly have independent phases, so that they fully 
account for the two unremovable phases in the model: Ninv,A = Nia = 2. The phase constraint 
from P22 33 ^^ t^^ same as in models 1 and 2 [^ The phase constraint from (5i3,32 is that the 
set of elements 

-^013,32 = {^13,^31,512,532} (5.3.4) 

cannot be made simultaneously real. The elements A31, ^33, and i?2i do not occur in any 
complex invariants and may be rephased freely. An example of an allowed form for the 
Yukawa matrices after rephasing of fermion fields is y'-"^ real. 



y(«)' 




(5.3.5) 



and 



yW 



/ 

I -821 1 

V 



\Bi2\e 
\B22\e 



A arg(Qi3^32) 
i argCPj'j' ) 



\ 



15 



32 I 





1^23 
1^33 I / 



(5.3.6) 



By our theorem on invariants and its specific Corollary 1, any 6'th order complex invariants 
have arguments which can be expressed in terms of those of the quartic complex invariants. 
We find that the model has one 6'th order complex invariant (so Ni^vfi = 1): 



^{ddu) 
V22,33,13 



i?22 533^135235^2^33 



(5.3.7) 



^^Ref. ||lj| actually presents a form in which the submatrix (5.1.7) of Y^'^'> is real. Such a form is not, in 
general, possible to achieve by any fermion rephasings, as we have already discussed above in the context of 
model 1. 
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whose argument is expressed in terms of those of the quartic invariants as 



argiQitslis) = arg(Pi2,33) + arg(Qi3,32) 



(5.3.8) 



The total number of complex invariants is Ninv = Ninv,A + Ni 



inv,6 



5.4 Model 4 



The fourth model is defined by 



y(") 



/ Ai2 

A21 A22 A23 

V A32 A33 



(5.4.1) 



/ Bi2 \ 

yW = B21 B22 (5.4.2) 

V 533 ; 

This model can be obtained from model 1 by the interchange of Y^'^^ and Y^'^K Using this 
fact, one can immediately determine all of the phase properties and invariants of the model 



from those of model 1: A^„ 



N. 



2, and the two independent quartic complex invariants 



are 



22,33 



^22^33^32^32 



(5.4.3) 



and Qi2,22 as given in ( [4.18 ). These fully account for all phase constraints in the model, 
which are given by those for model 1 with the interchange of F*^"^ '^^ Y^^'^ {^jk ^~^ -SjA:)- An 
allowed form for the Yukawa matrices after rephasing is 



y(")' 



and 



V 





I ^21 1 




\Ai2\e- 
\A22\e 



-i arg(Qi2,22) 
i arg(P2'2,33) 



\A 



32 



\ 

1^3 I 
1^331 / 



(5.4.4) 




(5.4.5) 



The model has one 6'th order complex invariant, <532'23 12 (given explicitly in ( ^.2.15|) for 
model 2). The argument of this 6'th order invariant is expressed in terms of those of the two 
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complex quartic invariants by the relation 

arg((33r23,i2) = - arg((5i2,22) - ^rgiP^^U) (5.4.6) 

As this illustrates, it is obvious that if the same 6'th order complex invariant occurs in 
two different model and if in both cases its argument is expressible in terms of those of 
the respective complex quartic invariants for these models, then these expressions will be 
different, if, as will generally be the case, the complex quartic invariants are different for the 
two models. Specifically, one sees that the relation ( p. 4. 61 ) in the present model differs from 
( p. 1.161) in model 2. 

The special case of model 4 in which \Y I \ = 1^^, |, f = u,d and F22 = ^22 = |l^ is 



not included in ref. [0 among its experimentally acceptable forms because of the large values 
of I V^fcl which it yields and, if other quantities are used as inputs to get a prediction for the top 
quark mass rrit, the high values of rrit which it predicts. However, this special case, and also 
its generalization with \Y ^ \ not necessarily equal to \Yy \ has several interesting properties 
concerning phases and rephasing invariants. First, since Neq is reduced to 9 while rank{T) 
remains at 8, there is only one unremovable phase, Np = 1. As is evident from eqs. ( |5.4.3[ ) 
and (f4.18|) , both of the quartic complex invariants P22 33 and Q12.22 vanish, so that there are 



no quartic complex invariants; Ninv,4 = 0. Therefore, in contrast to the general model 4 with 
nonzero ^22, in this special case, (i) the only complex invariant is the 6'th order invariant 
(532"23,i2 given above in ( p.2.15| ); (ii) as is evident from ( p.2.17| ), the argument of this 6'th 
order invariant cannot be expressed in terms of any lower complex invariants, since there 
are none. Thus this special case provides an example of a model in which the lowest-order 
complex invariant (s) are of 6'th rather than 4'rth order, and Ninv = Ninvfi = 1- The phase 
constraint yielded by this invariant is that the set of elements {^12,^23,^32,^33,512,-822} 
cannot be made simultaneously This allows y'-"^ to be made real; if it is, then {-612,-822} 
cannot be made real. The constraint also allows y*^"'^ to be made real, and if it is, then 
{y4i2, A23, A32, A33} cannot all be made real. 
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5.5 Model 5 



This model is given by 



y(") 



and F^'^) as in (^X^), 



yW 









^13 





^22 


^23 


^31 


^32 


^33 





Bl2 





B21 


B22 











-B33 



(5.5.1) 



Here N, 



eg 



10, and the T matrix is 10 x 9 with rank 8, so that Np 

(u) 



independent complex quartic invariants: -P22,33 ^^ ^^ (|5.4.3| ), and 



Ql3,22 — ^13-622^23-^12 



(5.5.2) 
2. We find two 

(5.5.3) 



Since the two complex quartic invariants have independent phases, they fully account for all 
of the phase constraints. These constraints are that it is not possible to render the respective 
elements comprising the sets 



Sr,(u 



(u) 



{A22, A23, A32, ^433} 



(5.5.4) 



or 



5'qi3,22 = {^13, ^23, -B12, -B22} (5.5.5) 

simultaneously real by any fermion rephasings. These constraints allow one to render Y^'^^ 
real. If one does this, then one of the two unremovable phases must reside in the submatrix 



A22 A23 
^32 ^33 



(5.5.6) 



formed by the elements of set ( |5.5.4| ) and the other in the set {A13, ^23}. The elements ^31, 
i?2i, and i?33 do not occur in any complex invariants and hence may be rephased freely. A 
possible form for the Yukawa matrices after rephasing is 



y{«)' 



/ 








14^2 g* arg(Qi3,22) \ 







/7-,(li) ^ 

A22 e' '''^'^^22,33) 


1^231 


\ 


1^31 1 


IA32I 


1^331 ; 



(5.5.7) 
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yW 



1 ° 


\Bl2 





B21 


-S22 





\ 





|533 



(5.5J 



Since the quartic complex invariants already contain all of the unremovable phases, it follows 
from our theorem on invariants and its Corollary 1 that any 6'th order complex invariant 
yields no new phase constraint. We find that there is one 6'th order complex invariant, 

^32,13,22 = ^32^13-822^33^22-^12 (5.5.9) 

Its argument is expressed in terms of those of quartic complex invariants as 

arg(QS,22) = - arg(Pi2,33) + arg((5i2,22) (5.5.10) 

The model thus has Ninv,4 = 2, Ninv,6 = 1, and Ninv = 3. 

5.6 Model 6 

We have noted above that, to the extent that one understands the phenomeno logical impli- 
cations of string theories, they do not appear to yield symmetric Yukawa matrices. This is, 
indeed, part of the reason that we have listed the previous five forms for these matrices in a 
form which is not, in general, symmetric. Here, as a further illustration of our methods, we 
determine the unremovable phases and associated complex invariants for a generalization of 
one specific model of Yukawa matrices inspired by a particular 4D superstring construction 
and sketched in P The model is defined by 



yW 



/ 


^12 





^21 


^22 


^23 


V 


^32 


^33 


/ 


B12 





521 


B22 


B23 


I 


B32 


B33 



(5.6.1) 



yW = B21 B22 B23 (5.6.2) 

\ B32 B33 J 

As can be seen, if one takes ^22 = 0, this reduces to model 4. The present model has 

Ngq = 12. We calculate the 12 x 9 T matrix as in the previous cases and find that it has 

^^Ref. M does not give any discussion of the number of unremovable phases or which elements of the 
Yukawa matrices can be made real. 
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rank{T) = 8, so that Np = 4. We find that there are Ninv,A = 7 quartic complex invariants, 
two of P4 type and five of Q4 type. These are ^22,33; -^22,33; Qi2,22, ^12,32, ^22,32, <522,33, 
Q22,33) and (523,33- The ordering of the arguments of the complex quartic invariants in the 
vector ^4 may be taken to be the same as in the list just given. The corresponding Z4 matrix 



for these order 4 invariants is the Nj 



invA 



X N, 



eq 



{Z^ 



4 6 



/ 1 



1 -1 
1 
1 
1 

\ 



-1 -1 












1 



7-row by 12-column matrix 












-1 
-10 

0-1 









1 
1 








-1 









\ 

-1 1 





1 
1 









1 



(5.6.3) 



(where the second subscript refers to model 6). We find that rank{{Z4)Q) = 4, so that the 
quartic complex invariants fully account for all of the Np = Nia = 4 unremovable phases. 
Removing the last three rows of (^4)5 does not reduce the rank, which shows that the 
arguments of the first four quartic complex invariants in (.^4)6 are linearly independent, so 
that we may take these complex invariants, ^22335 -^22335 Qi2,22, and (5i2,32j as the set with 
independent arguments. Each of these complex invariants implies a phase constraint that 
the corresponding set of elements cannot be made simultaneously real. A possible form of 
the Yukawa matrices, after rephasing, which is consistent with the phase constraints, is 



y(«)' 



/ 

V 





I ^21 1 




1^12 |e* ^''SWl2,22) 

1^22! 

I 4„„U«(arg(Qi2,22)-arg(Qi2,32)) 



I A 



23|e 



yW 



/ 

1-^21 1 

V 



I -^12 1 
I-B22I 
I-B32I 



\ 

j(-arg(P2'2J;3)-arg{Qi2,22)+arg(Qi2,32)) 
IA33I 



/ 



\B. 



23|e 



-i arg{P2*2,33) 



IB 



33 1 



(5.6.4) 
(5.6.5) 



/ 



The model has A^^ 



invfi 



4 independent 6'th order complex invariants: 



V22,33,12 



^{uud) 
V33,12,22 



^22^33-812^23^^2-^32 



^33^12-822^32^23-8^2 



(5.6.6) 

(5.6.7) 
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N (ddu) 



^22,33,12 as given explicitly before in ( |5.2.16|) for model 2, and 



^33,12,22 — -833^12^122-832-823^^2 (5.6.8) 

Thus, Ninv = Ninv,i + Ninvfi = 11 for this model. By our theorem on invariants and its Corol- 
lary 1, the arguments of all of these 6'th order complex invariants are expressible in terms 
of those of the four quartic complex invariants with independent arguments. The methods 
have been fully illustrated in our previous models, so we omit the explicit expressions. 



5.7 Model 7 

We proceed to mention two toy models briefly to illustrate certain theoretical points. We 
have exhibited a special case of model 4 for which the only complex invariant occurs at 6'th 
order and is of Qe type. It is of interest to exhibit a form for the Yukawa matrices which 
again has only a 6'th order complex invariant, but for which this is of Pq rather than Qe 
type. We have constructed such a model: 



y{") 



and r^'^) as in (|X|), 



yW 






A12 





A21 














^33 


5n 


-812 





-821 





-823 





-832 


-833 



(5.7.1) 



(5.7.2) 



This model thus has Nf.q = 9. We calculate the T matrix and find that it has rank{T) = 8, 
so that A'p = 1. By design, this model has no quartic complex invariants, and a 6'th order 



complex invariant 

-Pll, 23,32 = -811-823-832-812-821-833 (5.7.3) 

The graphical representation of this invariant is given in Fig. 3(b) for f = d. The elements 
of y^"-* are unconstrained by any phase invariants and, in particular, Y^'^^ can be made 
real. The phase constraint from -P112332 forbids Y'^'^^ from being made real or hermitian. A 
possible form for Y^'^^ after rephasing is 



yW 



/ l^iil 
I -821 1 

V 



I -812 1 


I -832 1 



\ 



I -823 1 e 



i arg{Pjj^ 23 32) 



\B: 



(5.7.4) 



33 



50 



5.8 Model 8 

We noted in section 2 that the rank of the T matrix is 8 for all of the realistic models that 
we have studied, but that, in principle, rank{T) may be less than 8. We present here a toy 
model to illustrate this. The model also provides an illustration of the inequality ( |4.23| ), 
that the number of unremovable phases in the Yukawa matrices is > the number of phases 
in the CKM quark mixing matrix, V. Although formally the model has Ng = 3, it is of a 
degenerate type in which the third generation is decoupled from the first two in both the 
up and down quark sectors. Hence, in some ways it acts more like an Nq = 2 model. The 
Yukawa matrices are 



So that Npa = 10. The T matrix is 




(5.8.1) 



(5.8.2) 



"eq 



/100100000\ 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 
10 10 

Vooioooooiy 



(5.8.3) 



with rank{Ts) = 7, so that Np 



N, 



eq 



rank{T) = 3. The model has Ni, 



N, 



invA 



complex invariants, which are precisely those listed in eqs. ( [4.14| ) for f = u and d and in 
eqs. (|4.15| )- (|4.18| ) in section 4 for the Ng = 2 model with all nonzero arbitrary elements. 
The Nia = Np = 3 complex invariants among these with independent arguments may be 
taken to be -Pii;22) Ai,22) ^^d (5ii,22- Note that this model also is like an A^^ = 2 model in 
that there is no CP violation. If one chooses to consider only the mutually coupled part of 
the model, then the Yukawa matrices would reduce to the upper left 2x2 submatrices of 
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( ^.8.1|) and (|5.8.2|) , the model would reduce to a quasi Nq = 2 model, and T would become 



an 8 X 9 matrix with rank 5. Hence for the mutually coupled sector of the model, the rank 
of T would obey the relation that rank{T) = rank{T)max = 3A^g — 1 (whence rank{T) = 5 
for Ng = 2) analogous to the rank 8 nature of the true Nq = 3 models. When regarded as a 
degenerate Nq = 3 case, this model also constitutes an exception to the empirical rule that 
Np = {Np)max — Nz, since {Np)max = 10 for Ng = 3 and A^^ = 8. However, again, if one 
considered only the mutually coupled sector of the model, then it would be a quasi Nq = 2 
model, so that {Np)max = 3, A^^ = 0, and the empirical relation would be obeyed. Such 
subtle distinctions are not necessary for the realistic models studied here. 

6 Conclusions 

The goal of understanding quark masses and quark mixing remains one of the most important 
outstanding problems in particle physics. When one constructs a model of the quark Yukawa 
matrices, it is necessary to determine how many unremovable phases there are in y(")) and 
Y^'^\ and to determine where these phases may be assigned, and which elements may be 
rendered real by rephasings of fermion fields. We previously reported a general solution 
to this problem [Q] involving several theorems and a related analysis of complex rephasing- 
invariant products of elements of the Yukawa matrices. In the present paper we have given a 
detailed discussion of our methods and results. We have also presented new results on phases 
and complex invariants for the case of arbitrary Nq, the number of fermion generations, since 
this elucidates their theoretical properties and gives further insight into the case of physical 
interest, Nq = 3. Finally, we have given a detailed application of our methods to currently 
viable models. 

This research was supported in part by NSF grant PHY-93-09888. We thank L. Lavoura 
for a discussion on the quartic P invariant. 
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